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An intriguing and elusive feature of the time-like hadron form factor is the possible presence 
of an imaginary part associated to rescattering processes. We find evidence of that in the recent 
and precise data on the proton time-like form factor measured by the BABAR collaboration. By 
plotting these data as a function of the 3-momentum of the relative motion of the final proton and 
antiproton, a systematic sinusoidal modulation is highlighted in the near-threshold region. Our 
analysis attributes this pattern to rescattering processes at a relative distance of 0.7-1.5 fm between 
the centers of the forming hadrons. This distance implies a large fraction of inelastic processes in 
pp interactions, and a large imaginary part in the related e + e“ —> pp reaction because of unitarity. 

PACS numbers: 


Electromagnetic hadron form factors (FFs) are funda¬ 
mental quantities which describe the internal structure 
of the hadron (for a recent review see Ref. Q)- FFs 
enter explicitly in the coupling of a virtual photon with 
the hadron electromagnetic current, and can be directly 
compared to hadron models which describe dynamical 
properties of hadrons. They are experimentally accessi¬ 
ble through the knowledge of the differential cross section 
and the polarization observables. 

Traditionally most data on FFs come from electron- 
proton elastic scattering. It is assumed that the inter¬ 
action occurs through the exchange of a single virtual 
photon which carries a four momentum transfer squared 
q 2 . In this kinematical region (space-like, SL) the virtual 
photon-proton coupling and the electron-proton scatter¬ 
ing cross section are described via two real FFs, electric, 
Ge, and magnetic, Gm- 

FFs have been also studied in the time-like (TL) re¬ 
gion of momentum transfer squared. They are measured 
through the reactions 

e+ + e _ -A p +p , (1) 

p + p -A e + + e~ , (2) 

where a hadron pair is formed by or annihilated into a 
virtual photon. In the following we will refer to the for¬ 
mer reaction, when not otherwise specified. 

Assuming single photon exchange, the unpolarized 
cross section contains the squared moduli of two TL FFs 
(electric and magnetic FF), which are complex-valued 
functions of q 2 . The imaginary part is expected to be 
large and information on the relative phase between Ge 
and Gm can be extracted from single spin polarization 
experiments !2] which are presently out of reach. In this 
letter evidence for periodic structures in TL FF data is 
reported and related to their complex nature. 

Many models for the hadron coupling to the virtual 


photon have been developed and applied to the calcula¬ 
tion of SL FFs. Some of them may be analytically con¬ 
tinued to the TL region. This is the case for approaches 
based on vector meson dominance @,3 and dispersion 
relations Constituent quark models in light front 

dynamics may be applied Q, as well as approaches based 
on AdS/QCD correspondence (0. A phenomenological 
picture has been recently proposed for an interpretation 
of FFs in both the SL and TL regions 0. 

The individual determination of the electric and mag¬ 
netic TLFF is obtained in principle, from the angular dis¬ 
tribution of reactions ClED, but until now the luminosity 
was not sufficient. The various experimental results are 
therefore compared on the basis of a generalized FFfl^. 
which is related to the unpolarized cross section by: er: 


\F P \ 2 


3/3q 2 a 


2ira 2 



(3) 


where a = e 2 /(4n), (3 = ^/1 — 1/r, r = q 2 / (4M 2 ), and 
M is the proton mass. 

Even in these simplified terms, it has long been diffi¬ 
cult to analyze with precision the behavior of the data 
over a broad kinematic range because of the uncertainties 
and of the matching of data from different experiments 
which covered limited (f_ regions. The recent data by the 
BABAR collaboration [111.Il2j| cover with reasonable con¬ 
tinuity a region ranging from slightly over the pp thresh¬ 
old to q 2 « 36 GeV 2 . In particular about 30 data points 
have been extracted in the region q 2 < 10 GeV 2 , with a 
relative error lower than 10 %. These features allow for a 
refined analysis of the systematic behavior of the TLFF, 
where large-scale and small-scale (in q 2 sense) properties 
of the data distribution may be scrutinized. 

From now on, we use the expression ” near-threshold re¬ 
gion” to indicate a q 2 -range extending from the threshold 
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of the pp channel up to q 2 ~ 9 GeV 2 (with the conven¬ 
tion c = h = 1). In this kinematic region, two different 
scales participate: the total energy of the colliding e + e~ 
pair is > 2 M « 1.9 GeV, while the kinetic energy of the 
created pp pair is relatively small. Therefore one may ex¬ 
pect to observe complex effects where a highly relativistic 
formation picture expressed in terms of quarks and glu¬ 
ons coexists with non relativistic interactions of two slow 
hadrons leaving the formation zone. 

Proton-antiproton interactions in the near-threshold 
region have been studied in experiments at LEAR (see 
[Hill and references therein for previous data) and 
more recently at AD [15]. These measurements could not 
separate spin channels and, as in the case of the single 
effective FF of Eq. ©, these data have been mostly an¬ 
alyzed in terms of a single effective scattering amplitude, 
as if proton and antiproton were scalar particles. 

We define ’’large inelasticity” when, writing the am¬ 
plitude as a sum of partial waves, at least half of the 
incoming flux is absorbed for all the partial waves of 
angular momentum L satisfying L < Rp , R ~ 1 fm. 
The unitarity limit is reached when there is total ab¬ 
sorption for these waves. For the inelastic cross section 
p + p —> X ^ p + p" large inelasticity” occurs in all the 
kinematical range of interest here, with the possible ex¬ 
ception of the region p < 50 MeV [3 Ell- For 100 
MeV the inelastic cross section of ~ 40 mb is close to the 
black disk limit ss 50 mb. 

As a consequence, unitarity leads to a large imaginary 
part in the amplitude of p + p —> exclusive final states, 
including ([2]). Near the threshold this is rigorously stated 
by Watson’s final state theorem [3 applied to reaction 
©• More in general the presence of a large transition 
amplitude p + p —> A n , where A n is an on-shell channel, 
implies a contribution to the imaginary part of the ampli¬ 
tude of 0 from a Cutkosky cut applied to the interme¬ 
diate state of the 2-step process p + p —>• A n —» e + + e~. 

Understanding at which extent phenomenological pp 
interactions could be final state interactions of the re¬ 
action © and seriously affect its amplitude, is however 
complicated by an evident ” mismatch” of two known fea¬ 
tures of these processes: 

1) The analysis of annihilation and scattering data has 
shown that colliding proton and antiproton do not over¬ 
lap at small kinetic energies. As soon as these particles 
come close within 1 fm of each other, they either annihi¬ 
late or scatter in a hard-core way [IQ, [l9|. So the wave 
function describing a pp state presents a ’’hole” of size 
1 fm. Within a momentum range of 400 MeV over the 
pp threshold, this property is demonstrated by counter¬ 
intuitive phenomena as the equality of pp, pD , and p 4 He 


annihilation cross sections at small energies and the 
suppressed effect of the electric charge in p -nucleus anni¬ 
hilation with the paradoxical effect of n cross sections on 
heavy nuclei exceeding p ones [§c| . 

2) In the exclusive reaction 0) quark-counting rules 


2ll [22] predict that in the initial stages of the forma¬ 
tion process quarks and antiquarks are concentrated in 
a region of size l/^/q 2 , which means at a relative dis¬ 
tance of 0.1 fm when q 2 ss 4 GeV 2 . So, in the initial 
formation stages the hadrons lie at a distance that is 
not normally tested in phenomenological pp interactions, 
making rescattering effects largely unpredictable. 

In order to search for signals of final state effects at 
small kinetic energies in the data, it is more convenient to 
introduce variables directly related to the relative motion 
of the hadron pair. In the following we will use the 3- 
momentum p of one of the two hadrons in the frame 
where the other one is at rest: 


= \JE 2 - M 2 , E = q 2 /(2M) - M. 


(4) 


The usefulness of this variable presumes that the process 
can be divided into two stages: formation and rescatter¬ 
ing, where the latter involves energies on a smaller scale 
than the former. This means that the amplitude for the 
process is the sum of a leading term due to a ’’bare for¬ 
mation” process taking place on a time scale 1/^fq 2 , and 
a relatively small perturbation associated to rescattering 
processes taking place on a larger time scale. 

A consequence of this assumption is that the measured 
FFs can be fitted by a function of the form : 

F(p) = F 0 (p) + F osc (p), (5) 

\Fosc(p)\ « |*o(p)|, (6) 

< F osc (p) > A p ->■ 0 for Ap » 1 GeV 2 , (7) 


where: 


• Fq ( p ) is the translation in terms of the variable p of 
a known parametrization that has been adjusted on 
the data in the full range 4M 2 < q 2 < 36 GeV 2 (see 
below) ignoring small-scale oscillations. Fq (p) is 
regular and smooth on the GeV/c scale. We name 
it “regular background fit”. 


• F osc (p) reproduces GeV-scale or sub-GeV-scale ir¬ 
regularities in the lower part of the p range. We 
name it “oscillation fit”. < F osc (p ) >a p is the lo¬ 
cal average of F osc (p) over the momentum range 
[p - Ap/2,p + Ap/2). 


The data by the BABAR collaboration 11, 12] are se¬ 
lected for this study, since they are the most precise data 
in the near-threshold region and they cover with conti¬ 
nuity a very large kinematic range. Both properties are 
necessary for our analysis. The fitting procedure is the 
one of the Minuit package of root.cern.ch based on 
the minimization of y 2 = Yli(f( x i>Pi) ~ 2/?:) 2 /of, where 
cji is the error on the point of coordinates (, Xi,yi ) and / 
is the fitting function depending on the parameters pt. 
The error on pi is the interval where y 2 /n.d./. increases 
by one unit ( n.d.f. is the number of degrees of freedom: 
number of points minus number of parameters). 
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Fig. 1: World data on the TL proton generalized FF as a 
function of q 2 , together with the calculation from Eq. m 
(blue dash-dotted line), Eq. © (red dashed line), Eq. (1101) 
(yellow long-dashed line), and Eq. 1121 (black solid line). The 
world data are from Ref. [Ill, 02] (black solid circles), Ref. 
[H (red stars), Ref. [24j (green squares), Ref. [27j| _ (blue 
triangles up), Ref. D/H ] (yellow triangles down ), Ref. [29|, [3Cj | 
(cyan full crosses), Ref. jR] (magenta full diamonds), and 
Ref. [32](dark green asterisk). The insert magnifies the near 
threshold region. 


To satisfy Eq. 0 a good choice of the regular back¬ 
ground term Fq( p) is needed. The generalized FF has 
been consistently extracted at e ± e" colliders and an¬ 
tiproton facilities. It shows a decreasing behavior as func¬ 
tion of q 2 , which was generally fitted in the experimental 
papers before the year 2006 with the function 24|, [25]: 


\Fscaling | — 


A 


(q 2 ) 2 log 2 (< 7 2 /A 2 ) 


( 8 ) 


A good fit of the data prior to BABAR last results was 
obtained with A = 40 GeV" 4 and A = .45 GeV 2 . 

Based on Ref. [H], in order to avoid ghost poles in a s 
the following modification was suggested (0): 


\Fscaling-\-corr{fl )| 


A 


(q 2 ) 2 [logV/A 2 )+7T 2 ] 


(9) 


Tab. I: Fit parameters from Eq. in. 


A ± AA 

B±AB 

[GeV]- 1 

C± AC 
[GeV] -1 

D±AD x 2 /n.d.f 

0.05 ±0.01 

0.7 ±0.2 

5.5 ±0.2 

0.03 ±0.3 1.2 


The TLFF data from the BABAR collaboration 11 


12] were obtained from the reaction 


P+P + 7) 


( 11 ) 


where the photon is preferentially emitted in the entrance 
channel. These data, extending from the threshold to 
q 2 ~ 36 GeV 2 , are well reproduced by the function 35] : 


\FBABAn(q 2 )\ 


_ A _ 

(1 + q 2 /m 2 ) [1 - g 2 /0.71] 2 


A = 7.7 GeV" 4 , m 2 = 14.8 GeV 2 . (12) 

The world data are shown in Fig. [T] as a function of 
the transferred momentum q 2 , and compared with the 
fits from Eqs. MUD - The near-threshold region is 
highlighted in the insert. 

In the following, we present results with Fo(p) = 
FBABAR[q 2 (p)\- This function does not follow the ex¬ 
pected asymptotic QCD counting rules, but best repro¬ 
duces the BABAR data, the slope of which is steeper than 
l/(q 2 ) 2 . It has to be considered as a local approximation 
of some more complicated function. We have checked 
that, taking any of the above background possibilities 
gives consistent results (within the errors) although the 
fit has a smaller y 2 /n.d.f. using Fbabar- 

In Fig. ©i the BABAR data are plotted as a func¬ 
tion of p. The result of the fit using Eq. 03 is then 
subtracted from the data. This difference V (i.e., data 
minus Fq( p)) is shown in Fig. [5b and exhibits a damped 
oscillatory behavior, with regularly spaced maxima and 
minima. Assuming that the first maximum is at p = 0, 
the distance between this, the 2nd and the 3rd maximum 
is A p Ri 1.14 GeV. After the 3rd maximum the oscilla¬ 
tions of the data are within the error bars. 

This behavior is fitted with the 4-parameter function 


F osc {p) = A exp {-Bp) cos (Cp + D). (13) 


In this case the best fit parameters are A = 72 GeV" 4 
and A = 0.52 GeV 2 . 

In Ref. [|J a form was suggested with two poles of 
dynamical origin (induced by a dressed electromagnetic 
current) 


|At3(<7 2 )| 


A 

(l-q 2 /m 2 )(2-q 2 /m 2 )- 


( 10 ) 


The best fit parameters are A = 1.56, m 2 = 1.5 GeV 2 
and m 2 = 0.77 GeV 2 . 


The values of the parameters are reported in Table [H 
The relative size of the oscillating term over the reg¬ 
ular background is ~ 10 %. The damping range of the 
oscillations of Fig. [5b is 1/B « 1.4 GeV. F 0 (p) decreases 
by a factor 1/e within about 1.5 GeV. The relative mag¬ 
nitude of the oscillations to the regular background term 
Fq(p) does not change much at increasing p , although 
increasing errors make the oscillations undetectable for 
p > 3 GeV (q 2 > 10 GeV 2 ). At asymptotically large 
q 2 values, the Phragmen-Lindeloff theorem [0 requires 
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Fig. 2: (a): TL proton generalized FF as a function of p from 
Ref. 11, []J] ; the line is the regular background fit with Eq. 
(11211 : (b): data after subtraction of the fit; the line is the fit 
with Eq. m- 


that the imaginary part of TLFFs vanishes, which implies 
that rescattering disappears. So, although we expect a 
large-momentum suppression of the relative weight of the 
rescattering terms, this is not seen in the range q 2 < 10 
GeV 2 where error bars allow us to distinguish systematic 
from statistical oscillations in the data. 

The periodicity and the simple shape of the oscillations 
seem to exclude a random arrangement of maxima and 
minima of heterogeneous origin. Rather, they indicate 
a unique interference mechanism behind all the visible 
modulation. A simple oscillatory behavior in p means 
that the waves corresponding to the outgoing particles 
originate from a small number of coherent interfering 
sources: these waves may share a common initial source 
but be rescattered along different paths or in different 
ways so to acquire different phases. We may speak of 
“alternative rescattering pathways”. These must be in 
a small number and there must be discontinuity among 
them, otherwise we would see diffraction patterns instead 
of interference patterns. For example, in Q it has been 
suggested that, during the intermediate stages of pp for¬ 
mation, charge and color are distributed in a highly in¬ 
homogeneous way, with discontinuity features. 

Since we do not know the rescattering mechanism, we 
cannot identify the sources of rescattered waves, but we 
may gain some clue on their space distribution. Let f 
be the space variable that is conjugated to p via three- 
dimensional Fourier transform. We may identify r as 
the distance between the centers of the two forming or 
formed hadrons, in the frame where one is at rest. Let 
Mo(r) and M(r ) be the Fourier transforms of the regular 



Fig. 3: (Left) JWo(r), Eq. (fl4l) . (Right): M(r) (solid line) 
from Eq. m, and Mo(r) (dashed line) for comparison (linear 
vertical scale). 


background fit and of the complete fit: 

Fo(p) = J d 3 r exp(ip ■ r) M 0 (r) (14) 

F{p) = F 0 (p) + F oac (p) = J d 3 r exp(ip ■ r)M(r). (15) 

Mo(r) is shown in Fig. [3] left panel. The most relevant 
feature is that M 0 (r) decreases by 7 orders of magnitude 
for r ranging from 0 and to 2 fm. The decrease is reg¬ 
ular and almost constant on a semilog scale. Mo(r) is 
steep near the origin, too. From a mathematical point 
of view this follows from the fact that at the threshold 
of the pp channel the function Fq(p ) is a regular, contin¬ 
uous, and rapidly decreasing function of q 2 . It can be 
interpreted by the fact that both Fq (p) and its transform 
M 0 (r) are exp ression of that short distance quark-level 
dynamics [ 2 lL [ 221 ] that permits exclusive pp production 
at the condition that the final quarks and antiquarks are 
formed within a small region. Near threshold, the size of 
this region is < 0.1 fm, much smaller than the standard 
hadron size. 

In the right panel of Fig. [3 \M{r) is superimposed to 
M 0 (r). We notice that these two functions do not dif¬ 
fer for r < 0.7 fm, and that the physical reason of the 
data oscillation must be searched for in processes taking 
place in the r-range 0.7-1.5 fm. This range is important 
because it includes the distances corresponding to the 
largest annihilation probability in the phenomenolo gica l 
pp interactions in the near-threshold region [SEUi 
At a distance of 1 fm, the relevant part of rescattering 
must involve physical or almost physical hadrons that an¬ 
nihilate into groups of 2-10 mesons. As discussed above, 
this means a large contribution to the imaginary part 
of the amplitude for 7 * —>• pp from the Cutkosky cuts 
applied to all the 2 -step processes like 7 * — mr —> pp, 
where mr is a state composed by n on-shell pions (or 
other mesons). 
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Rescattering with a phase shift between alternative 
channels may take place via formation of s-channel poles, 
or via t —channel photon/meson exchanges. Both classes 
of processes have important SL continuations. Phe¬ 
nomenological s —channel poles lead to nontrivial phe¬ 
nomena due to their imaginary parts, but not to periodic 
oscillations (see the analysis in 37] on the continuation 
of some known SL models to the TL region). Exchange 
of i-channel photons leads to a Coulomb phase in the TL 
region. In the case of photon/meson f-channel exchange, 
we have a large set of possibilities (e.g. nn formation 
followed by charge exchange). In the corresponding SL 
diagrams a virtual photon/meson is emitted by the nu¬ 
cleon before the hard vertex and reabsorbed after it. In 
the case of virtual photons, at large (TL and SL) |< 7 2 | this 
process would not be more relevant than multiple pho¬ 
ton exchange between the nucleon and the lepton cur¬ 
rents [38j. If the exchanged bosons are mesons, the cor¬ 
responding SL diagrams modify the distribution of the 
proton charge. The way SLFFs are affected is strongly 
model dependent as shown in particular in the neutron 
case in a series of works (see j39( and references therein). 

In summary, a systematic modulation pattern in the 
TLFF measured by the BABAR collaboration in the 
near-threshold region has been highlighted in the range 
q 2 < 10 GeV 2 . This modulation presents periodical fea¬ 
tures with respect to the momentum p associated with 
the relative motion of the final hadrons. It suggests an in¬ 
terference effect involving rescattering processes at mod¬ 
erate kinetic energies of the outgoing hadrons. Such pro¬ 
cesses take place when the centers of mass of the pro¬ 
duced hadrons are separated by ~ 1 fm. For this rea¬ 
son at least a relevant part of rescattering must con¬ 
sist of interactions between phenomenological or almost 
phenomenological protons and antiprotons. These phe¬ 
nomenological reactions are known to have inelastic cross 
sections overcoming 1/2 of their unitarity limit. Unitar- 
ity arguments imply the presence of a large imaginary 
part of TLFF. The relative errors of the data increase 
with q 2 , making us able to detect the modulation for q 2 
< 10 GeV 2 , but its relative magnitude of about 10 % is 
constant in this range, suggesting the interesting possi¬ 
bility that this modulation could be observed at larger 
q 2 in forthcoming more precise data. Precise measure¬ 
ments in the near threshold region are ongoing at BE- 
SIII (BEPCII), on the proton as well as on the neutron, 
bringing a new piece of information. The measurement 
of TL FFs in a large q 2 range will be possible at PANDA 
(FAIR). 
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